Enlarged Molecules from Excited Atoms in Nanochannels 
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The resonance interaction that takes place in planar nanochannels between pairs of excited state 
atoms is explored. We consider interactions in channels of silica, zinc oxide and gold. The nanosized 
channels induce a dramatically different interaction from that in free space. Illustrative calculations 
for two lithium and cesium atoms, demonstrate that there is a short range repulsion followed by 
long range attraction. The binding energy is strongest near the surfaces. The size of the enlarged 
molecule is biggest at the center of the cavity and increases with channel width. Since the interaction 
is generic, we predict that enlarged molecules are formed in porous structures, and that the molecule 
size depends on the size of the nanochannels. 
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How molecular interactions and energy transfer are af- 
fected by the shape and size of cavities and modulated 
by the surfaces that confine them is of interest in catal- 
ysis. Little work has been done in this area. Zeolites, 
mesoporous silicas, and their use for oil cracking may be 
one of the largest of all industries. Despite this, there is 
still no general agreement on the mechanisms of adsorp- 
tion of enzymes and other molecules in mesoporous silica, 
and their excited state interactions. An investigation 
of how molecular interactions and energy transfer are af- 
fected by the shape and size of cavities and of the surfaces 
that confine them is then of interest. Hopmeier et al. [2| 
have demonstrated experimental evidence for enhance- 
ment of dipole-dipole interaction in a microcavity and 
Agarwal and Gupta |3( have demonstrated this in theory. 
An inhibition to such work has been that the standard 
theoretical expression for the resonance interaction be- 
tween excited state-ground state atoms is incorrect [3-Q] ■ 

In this work we demonstrate how resonance interac- 
tions between excited atoms are strongly modified at 
nanoscale dimensions when the atoms interact inside pla- 
nar channels. We show that the containment effects 
on the interaction can lead to the formation of pecu- 
liar enlarged molecules. As compared to our previous 
work [j| the present contains a deeper analysis of the phe- 
nomenon. This includes an account for the origin of the 
short-range repulsive and long-range attractive interac- 
tion via spectral plots of interactions from different ex- 
citation branches and detailed studies of the effects due 
to different locations of the atomic species and different 
cavity sizes. The binding energy is strongest near the 
surfaces. The size of the enlarged molecule is biggest 
in the center of the cavity. We use lithium and cesium 



atoms and channels in silica, zinc oxide and gold as ex- 
amples. We first briefly rehearse the (correct) theory of 
the resonance interaction energy in channels and in free 
space. With that established we present some illustra- 
tive results. We compare the very different interactions 
of atoms in free space and in nanochannels. 

We have shown previously 0, [H] that, due to too dras- 
tic approximations, the underlying theory of resonance 
interactions derived from perturbative quantum electro- 
dynamics (QED) is only correct in the non-retarded limit. 
To see this we recall the standard argument: Consider 
two identical atoms where one initially is in its ground 
state and the other is in an excited state. This state 
can also be represented by a superposition of states: one 
symmetric and one antisymmetric with respect to inter- 
change of the atoms. While the symmetric state is likely 
to decay into two ground-state atoms, the antisymmet- 
ric state can be quite long-lived. The system can thus 
be trapped in the antisymmetric state 0, 3 ■ The energy 
migrates back and forth between the two atoms until ei- 
ther the two atoms move apart or a photon is emitted 
away from the system. First order dispersion interac- 
tions are caused by this coupling of the system, i.e. the 
energy difference between the two states is separation (p) 
dependent. After writing down the equations of motion 
for the excited system it is straightforward to derive the 
zero temperature Green function for two identical and 
isotropic atoms [1, [H, Q . The resonance frequencies (w r ) 
of the system are given by the following equation Q : 



1 - a{l\u)a(2\uj)T(p\uj) 2 = 0, 



(1) 



where the atom polarizability is taken to have the form 
(for lithium and cesium) 
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The parameters for Li and Cs were obtained [10] through 



2 




p (units of a ) p (units of a ) 



Figure 1. (Color online) Full interaction energy between two 
lithium atoms (when the modes in the (±;a;) branch are ex- 
cited) and various limiting results as functions of separation, 
p. The full result from Eq. @ is shown as the thick solid curve 
with circles; the large separation asymptote from Eq. <I10[) is 
represented by a thin curve with vertical bars; the T = K 
version of this asymptote from Eq. pip is represented by the 
steepest of the two dotted straight lines; the other dotted line 
is the n = term from Eq. (JT^J); the non-retarded full result 
is given by the solid straight line. All results are in atomic 
units, i.e., the energies in Hartree units and the separation in 
Bohr radii. 



a fit of the van der Waals and Casimir asymptotes to the 
results from a full ab initio quantum mechanical calcula- 
tion by Marinescu and You [11| . The fitting procedure for 
this single oscillator model works very well for all com- 
binations of alkali-metal dimers in vacuum. For dimers 
immersed in a medium one may expect that a more re- 
alistic form of the polarizabilities are needed. 

In a vacuum the field susceptibility matrix, T(p\iu>), 
has the following non-zero matrix elements [||, 



1 V 



(3) 



when the atoms are aligned along the x-axis. The sus- 
ceptibilities in a cavity with dielectric walls were used by 
Rahmani et al. 1 1 2| to discuss fluorescence lifetimes. If 
the two molecules are at the distance z from one of the 
channel boundaries the corrections to the free space re- 
sults when the atoms are symmetrically excited involve 
the change in the trace of the matrix, T^ x + T yy + T] z , 
where 
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Figure 2. (Color online) Full interaction energy between two 
symmetrically excited lithium atoms and between two sym- 
metrically excited cesium atoms at the center of a 1 nm planar 
channel in silica. 



Here Jo(qp) is a Bessel function of the first kind, 7, 



^fq 2 -\- 6iU) 2 /c 2 and 



A, 



7i - 7o£i 



A S = H-I^. (5) 
7i + 7o£i 7i + 7o 

The resonance energy in narrow channels can be cal- 
culated if the dielectric function, ei(iw), of the wall ma- 
terials is known. In the present work we need the di- 
electric functions for silica (SiC^), zinc oxide (ZnO), and 
gold ( Au) . The complex dielectric functions of the oxides 
were therefore determined employing a first-principles 
approach within the GW method based on the density 
functional theory [13j . The dielectric function of gold was 
determined from optical data [l4| . 

In the case of two identical atoms the above resonance 
condition can be separated in one antisymmetric and one 
symmetric part. Since the excited symmetric state has 
a much shorter life time than the antisymmetric state 
the system can be trapped in an excited antisymmetric 
state [J| . The resonance interaction energy of this anti- 
symmetric state is, 



U (p) = h[u> r (p) — cj r (00)] . 



(6) 
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Since the relevant solution of Eq. (TT]) really is the pole 
of the antisymmetric part of the underlying Green func- 
tion we can in a standard way [14] deform a contour of 
integration around this pole to obtain a both simple and 
exact expression for the resonance interaction energy, 

/>OC 

U(p) = (H/n) d^n[l + a(l\iOT(p\iO]. (7) 
Jo 

To account for the temperature (T) dependence we 
simply replace the integration over imaginary frequencies 
by a summation over discrete Matsubara frequencies [Til 
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f„ = 2irk B Tn/h, (8) 
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Figure 3. (Color online) Contribution from each term n for 
the system shown in Fig. [5] for a case when the total interac- 
tion energy is repulsive (p = 20ao). The frequency decomposi- 
tion of the interaction energy from T xx , T yy , T zz contributions 
and from the sum of these surface corrections, U xx , U yy , U zz 
and U^ at i — U xx + U yy + Ul z , respectively, is shown. Also 
shown is the total interaction energy contribution, t7 t(rfaI , in- 
cluding both surface corrections and the free space result from 
different frequencies. The discrete frequencies are shown at 
the top and at the bottom as circles. Note that for illustrative 
purposes the zero frequency contribution has been moved to 
the vertical axis. 



where ks is the Boltzmann constant and the prime in- 
dicates that the n — term should be divided by 2. 
Focusing first on what happens in free space we find that 
the correct leading term, at large separations, when the 
modes in the (±;x) branch are excited, is 

U(p,T) ~±^^jr / a(it n )e- xn {l + xn + x 2 n 2 }, (9) 

" 71=0 

where x = 2nkBTp/(hc). We can replace the polariz- 
ability with a(0) at intermediate and large separations. 
Within this approximation the resonance free energy is, 

x [1 + e 3x - e x (l + 2x - 2x 2 ) + e 2x (-l + 2x + 2x 2 )]. 

For intermediate values of x (see the middle portion of 
the full result in Fig.Q]) this free energy of resonance in- 
teraction varies as l//? 4 : 



U(p) ~ ±4hca(0)/(irp 4 



(11) 



This represents the dominating term in the interaction 
energy for oscillators in an excited configuration at zero 
temperature. However, for any finite temperature at 
sufficiently large distances the long-range interaction is 
dominated by the n — term. This term is here 



U(p,T) n=0 = ±k B Ta(0)/p 3 



(12) 



Figure 4. (Color online) Same as Fig. [3] but for a larger sepa- 
ration (p = 40ao ) between the lithium atoms so that the total 
interaction energy is attractive. 



This is the correct asymptotic long-range resonance in- 
teraction at any finite temperature (see the rightmost 
portion of the full result in Fig.[T]) 

As a first illustration of the full result and of several 
limiting expressions we present, in Fig.[TJ the results for 
two lithium atoms in free space when the modes in the 
(±;x) branch are excited. 

An interesting case studied in Fig.[5]is when the lithium 
or cesium atoms are excited symmetrically (i.e. when 
they are averaged over all possible orientations into a 
symmetric final state). The total resonance interaction 
produces short range repulsion and long range attraction 
creating the possibility to have bound states. In other 
words very large molecules can be formed inside pores. 
As we show in the figure the attractive potential well is 
about twice as deep when lithium atoms are replaced by 
cesium atoms. 

We show in Figs. [3] and|4] how different modes con- 
tribute to a total repulsive short range interaction and 
a long range attraction. In the free space case attrac- 
tive and repulsive interactions from the different branches 
cancel out in the non-retarded case. However, a different 
scenario occurs when the modes are confined. The differ- 
ent branches are squeezed out differently when the atoms 
are within a planar cavity. The x branch is directed be- 
tween the atoms while the z branch is perpendicular to 
the plates (the y branch is perpendicular to the z and x 
branches). The short range repulsion between two atoms 
in the cavity comes from the repulsive y and z branches 
being larger than the attractive x branch. As the distance 
between the atoms increases only the y branch gives re- 
pulsion and the total interaction is attractive. 

The depth of the attractive potential well is much 
larger if the atoms move away from the center of the 
channel and approach one of the boundaries. This is il- 
lustrated in Fig.(S] The attractive potential well is much 
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Figure 5. (Color online) Full interaction energy between two 
symmetrically excited lithium atoms in a 1 nm planar channel 
in silica (dashed curves), in zinc oxide (solid curves) and in 
gold (dotted curves). Here the atoms are positioned at the 
distance z from one of the walls. 
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Figure 6. (Color online) Minus the full interaction energy 
between a pair of lithium atoms in a symmetrically excited 
final state between two gold surfaces distances d = f nm and 
d — 2nm apart. 

larger for atoms in between gold surfaces as compared 
to between silica or zinc oxide surfaces. In Fig.[S] we 
show that the binding energy increases as the atoms move 
closer to one of the gold surfaces. The size of the molecule 
is largest at the center of the cavity and increases with 
channel width. However, for too large cavities the bind- 
ing energy in the center will be too small to effectively 
bind the atoms together. Close to surfaces, independent 
of channel width, there will be sufficient attractive inter- 
action to bind atoms together. 

The focus of this paper is on resonance interactions 
when excited atom pairs are inside nanochannels, e.g. 



inside porous silica or between gold walls. Resonance in- 
teractions can be suppressed inside channels when the 
excitation is not symmetric. For a symmetric excitation 
in free space the resonance interaction energy averages to 
zero in the non-retarded limit. A very different scenario 
evolves inside a planar nanochannel. Here there can be 
short range repulsion followed by long range attraction 
creating a possibility to have a bound state. This means 
that there can be very large molecules with a size that 
depends on the channel width. The binding energy is 
dramatically enhanced close to one of the two surfaces. 
Casimir-Polder interaction between atoms and surfaces 
drives the pair towards the surface where the size of the 
molecules decreases. The analysis of the resonance inter- 
actions undertaken here and its application to our par- 
ticular case of atoms in a nanochannel show some new 
features of interactions in confined geometries. Interac- 
tions in, e.g., cylinders can be tackled straightforwardly 
by available semi classical techniques [15j . Already, it 
is clear from the example we have studied that confine- 
ment geometry at nanoscales can produce qualitatively 
new results. We have confined the study to two atoms 
at the same distance from the surfaces. An extension 
of the work would be to investigate how the effects are 
modified when one atom is closer to a surface than the 
other. One knows that the decay rate of a dipole changes 
in a cavity and depends on the orientation [3] . The for- 
mation of these large scale molecules requires two atoms 
where one is in an excited state. This could be difficult 
to achieve in a cavity with thick walls especially if they 
are made of gold. One way to manage this for a metallic 
cavity could be to excite the atoms outside, before they 
enter the cavity. Another way might be to have one wall 
made of glass coated by a thin ITO (indium tin oxide) 
layer or of ZnO where the inner surface layer is heavily 
doped; this wall would be transparent in the visible and 
ultra violet parts of the spectrum and behave as a metal 
for lower frequencies. Then the atoms could be excited 
with a laser through this prepared wall. To describe this 
an extended Green tensor for stratified media has to be 
used. This is available in the literature [l6| . 

The observed effects could possibly be studied and ver- 
ified using gas chromatography with new filter types. 

Finally we may speculate in some possible applica- 
tions of the observed phenomena: design of more effective 
lubricants; better selection of molecules for "electronic 
noses"; more effective surface reactions for H2 generation; 
pressure sensitive filtering of gas flow in nanomaterials; 
new materials for gas chromatography. 
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